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We synthesize sub-THz longitudinal quasi-monochromatic acoustic phonons in a SrTiOs single 
crystal using a SrRu03/SrTiC>3 superlattice as an optical-acoustic transducer. The generated acous- 
tic phonon spectrum is determined using ultrafast X-ray diffraction. The analysis of the generated 
phonon spectrum in the time domain reveals a k-vector dependent phonon lifetime. It is observed 
that even at sub-THz frequencies the phonon lifetime agrees with the 1/w power law known from 
Akhiezer's model for hyper sound attenuation. The observed shift of the synthesized spectrum to 
the higher q is discussed in the framework of non-linear effects appearing due to the high amplitude 
of the synthesized phonons. 
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I. INTRODUCTION 

The increasing importance of coherent phonon spec- 
troscopy in material science is related to the growing 
problem of heat dissipation in modern nanoscale de- 
vices. This problem is impossible to solve without de- 
tailed understanding of underlaying phonon-phonon and 
phonon-electron interactions on the nanoscale. One of 
the methods to study this processes is coherent phonon 
spectroscopy, in which a particular phonon spectrum is 
excited coherently in the sample and detected optically. 
Research efforts in this direction resulted in significant 
progress in generation and detection of coherent phonons 
in various materials. The available phonon frequency 
has reached the THz acoustic limit 1 and basically the 
whole phonon frequency range nowadays could be ex- 
cited coherently. However, convenient optical detection 
methods based on Raman 2,3 or Brillouin 4 scattering al- 
low for the observation of phonons excited only in the 
vicinity of the Brillouin zone center. Therefore, sub- 
THz acoustic phonons could be accessed optically only 
in multilayer structures, in which the acoustic disper- 
sion branch backfolds many times inside a mini-Brillouin 
zone of a multilayer 5 . Modern progress in pulsed laser 
techniques as well as in multilayer fabrication has lead 
to a set of successful experiments in which the coherent 
zone-folded superlattice phonons have been optically ex- 
cited and detected 5-7 . However, these optical methods 
are insensitive to the THz frequency phonons which have 
propagated into the bulk of the crystal due to the unfold- 
ing of the phonon dispersion curve. Convenient optical 
methods based on Brillouin scattering in this case have 
a detection limit in the 100 GHz range given by the wave 
vector magnitude of the optical light 8 ' 9 . Recently, ultra- 
fast X-ray diffraction (UXRD) has become available to 
extend the accessible phonon frequency range to above 
100 GHz. It has been used successfully to study both the 



time-domain structure of optically excited zone-folded 
coherent acoustic phonons in epitaxial multilayers 10 as 
well as to observe the propagation of unfolded phonons 
into the bulk 11 . 

In this paper we report our new UXRD experiments 
from coherent quasi-monochromatic longitudinal acous- 
tic phonons in SrTiOs synthesized by fs-laser excita- 
tion of SrRu0 3 /SrTi0 3 (SRO/STO) epitaxial multilay- 
ers. Using UXRD we determine the laser excited phonon 
spectrum in SrTiOs and monitor the modification of the 
spectrum in the time domain. The epitaxial multilayers 
were prepared using Pulsed Laser Deposition 12 . The ex- 
periments are carried out at the BESSY EDR bcamlinc 
using a unique setup for a 1 MHz repetition rate UXRD 
experiments. 

The experiments are done in a traditional scheme 
which uses an optical delay line to change the time inter- 
val between the optical pump and the X-ray probe pulses. 
We use infrared optical pulses with the wavelength of 
1.03 /im for pumping and 8keV X-rays for probing the 
lattice dynamics. The important feature of this setup is 
the simultaneous acquisition of the X-ray photons scat- 
tered from the sample before and after the pumping op- 
tical pulse. This makes the X-ray intensity difference 
signal sensitive only to those changes in the crystal lat- 
tice which were exclusively initiated by the optical pulses. 
For a further details we refer to a recent publication de- 
scribing the setup 13 . 



II. THEORY 

A. Synthesis of quasi-monochromatic coherent 
acoustic phonons 

Recent studies showed that the optical excitation 
of a metal transducer by a sequence of ultrashort 
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TABLE I. Structural and mechanical properties employed in 
the calculations. 
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FIG. 1. (Color online) (a) Sketch of the 5-period epitaxial 
SrRu03/SrTiC>3 multilayer on top of the SrTiOs substrate. 

(b) Calculated instant strain profile in the sample along the 
surface normal, taken at 250 ps after the multilayer excitation. 

(c) Calculated spectrum of the strain pulse in the wave vector 
and frequency domains induced by the optical excitation. The 
vertical scale stands for the excitation fluence of 6mJ/cm 2 
used in the experiment. 



laser pulses is an efficient method to generate sub- 
THz quasi-monochromatic longitudinal acoustic (LA) 
phonons 14-16 . In essence, the repetitive generation of 
bipolar strain pulses by the laser-excited transducer 17,18 
forms a phonon wavepacket of narrow spectral bandwidth 
propagating throughout the substrate. In this report we 
consider a different approach which uses a spatial rep- 
etition instead of a temporal one, i.e., the excitation of 
a periodic metal-dielectric multilayer (superlattice) with 
a single ultrashort laser pulse (see Fig. 1(a)). This way 
the so-called superlattice phonon mode is excited 10, 19-22 
which subsequently unfolds into the substrate thereby 
forming LA phonon wavepackets with similarly narrow 
spectral bandwidth 11 ' 21 . In most cases the laser-induced 
coherent lattice dynamics may be calculated using either 
a model of a continuous elastic medium 17 or a linear- 
chain model (LCM) of masses and springs 21 . Here we 
employ the latter approach which will be more appro- 
priate for short-period superlattices and automatically 
accounts for the acoustic phonon dispersion. The details 
of the numerical model can be found in Ref. 21. The 
excellent agreement of such a linear-chain model with re- 
lated optical and UXRD experiments has been previously 
demonstrated for various sample structures 16, 18 > 20 ~ 24 . 

Further we consider linear chain calculations for a 
five-period SRO/STO multilayer with the spatial pe- 
riod D = 140 nm which is schematically shown in 
Fig. 1(a). The structural parameters of the sample which 
are determined by static high-resolution x-ray diffraction 
(HRXRD) 25 are collected in Table I. In addition, the 
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longitudinal sound velocities of the individual materials 
are shown 26 . The results of the linear-chain calculations 
using the parameters given in Table I and the experimen- 
tal pump fluence of 6 mJ/cm 2 are shown in the Figs. 1(b) 
and (c). 

The graph in Fig. 1(b) shows the calculated one- 
dimensional strain profile in the sample 250 ps after the 
excitation. Figure 1(c) plots the spectral amplitude of 
the linear-chain eigenvectors (normal modes) as a func- 
tion of the eigenfrequency. This amplitude spectrum is 
solely determined by the initial conditions 21 . For the 
bulk STO substrate the eigenvectors are plane elastic 
waves with wave vector q satisfying the well-known dis- 
persion relation of acoustic phonons 27 . As Fig. 1(b) il- 
lustrates, the optical excitation of the metal layers of 
the multilayer system results in the generation of a co- 
herent strain wavepacket propagating into the STO sub- 
strate at the longitudinal sound velocity 21 . The result- 
ing wavepacket inside the STO substrate attains the 
particular shape shown in Fig. 1(b), namely, five lead- 
ing compression pulses and five trailing expansion pulses 
which are separated by A ~ 140 nm, respectively. In 
other words, the metal/dielectric multilayer acts as the 
photoacoustic transducer synthesizing the coherent LA 
phonon wavepacket in the STO substrate. The sharp 
static profile of the thermal strain inside the multilayer 
(see Fig. 1(b)) remains unchanged with time because the 
linear chain model neglects the effect of heat diffusion. 
The heat diffusion in multilayers is a complicated sep- 
arate topic which lays out of the scope of this article. 
In this article we focus on the coherent lattice dynamics 
in the STO substrate which occur at a later timescale 
when the effect of heat diffusion within the multilayer 
does not play a role. For the detailed description of the 
wavepacket strain profile and its generation we refer to 
our earlier works 21 . 

The calculated amplitude spectrum in Fig. 1(c) con- 
tains several equidistant peaks. The most pronounced 
peak at go = rad/nm is responsible for the overall bipo- 
lar shape of the wavepacket 17,18 ' 21 . The width of the peak 
is deterined by the total thickness of the multilayer (Aq ss 
7r/5A). The peak around qi = 2-7r/A « 0.046 rad/nm 
corresponds to the characteristic spatial period A of the 
wavepacket. The non-sinusoidal shape of the wavepacket 
gives rise to the higher harmonics at integer multiples of 

Altogether, we find that using a periodic metal- 
dielectric multilayer as photoacoustic transducer we can 
generate LA phonon wavepackets similar to the quasi- 
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monochromatic wavepackets produced by multiple-pulse 
excitation of a thin metal film 14-16 . In both cases 
the wavepackets exhibit narrow spectral bandwidth and 
higher harmonics of lower amplitude. 

B. Ultrafast X-ray diffraction from sub-THz elastic 

waves 

The X-ray diffraction from crystals which are subject 
to a strong acoustic field is a well-established topic 28-30 . 
However most of the previous studies deal with strain 
fields generated by surface acoustic wave (SAW) trans- 
ducers. Such devices normally generate acoustic waves 
with wavelengths longer than either the X-ray extinction 
length or the X-ray coherence length. The X-ray diffrac- 
tion from a crystal lattice perturbed by such waves re- 
sults in modifications of the Bragg peak shape within the 
Darwin width or in the appearance of diffuse scattering 
contributions in the vicinity of the peak 31 ' 32 . The de- 
scription of the X-ray scattering from such waves usually 
requires dynamical X-ray diffraction theory. 

Here we consider quasi-monochromatic coherent LA 
phonons which in fact are elastic waves at hypersonic 
frequencies. The corresponding wavevectors q are large 
enough to allow for coherent Bragg-like scattering of X- 
rays from the associated "moving gratings". Due to the 
sufficiently high q vectors of the quasi-monochromatic 
phonons the X-ray scattering contributes in the off-Bragg 
region in which the X-ray scattering efficiency from the 
bulk of the crystal is small. This allows to probe the X- 
rays exclusively scattered from the wavepacked with only 
minor perturbation by the bulk-scattered X-ray wave- 
field. 

In this section we introduce the necessary theoreti- 
cal basis which allows a thorough interpretation of the 
UXRD experiments. Since we intend to apply the fol- 
lowing theory to study the laser-induced structural dy- 
namics in one dimension we restrict ourselves to a one- 
dimensional formulation. 

A plain elastic wave with wave vector q contributes 
to the scattered X-ray intensity if the X-ray scattering 
vector Q is given by 

Q = G±q (I) 

where G is a reciprocal lattice vector and Q = |k - kj| is 
the X-ray scattering vector 16,33-35 . 

From simulations of the scattered X-ray intensity using 
dynamical theory of X-ray diffraction one finds that the 
X-ray intensity scattered from the crystal perturbed by 
a bunch of elastic waves can be well described by the 
equation 

(I p (Q))t = Iu P (Q) + aA(q) 2 (2) 

where I p (Q) and I up (Q) is the scattered X-ray inten- 
sity from the perturbed and unperturbed crystal, respec- 
tively. The angle brackets stand for time averaging. The 



function A(q) is the spectral amplitude of the elastic wave 
with wave vector magnitude q = \Q — G\ and a is some 
constant. It is worth to show here that formula (2) is 
equivalent to the expression describing thermal diffuse 
scattering (TDS) from acoustic phonons 36 . To show this 
we need to relate the energy of a classical plane elastic 
wave in the crystal with the phonon population. The en- 
ergy of plane elastic waves in the classical linear theory 
of elasticity is proportional to the squared product of the 
wave amplitude A and frequency u> 

E{uj) oc A 2 lo 2 (3) 

In a crystal lattice this corresponds to the energy of the 
corresponding vibrational normal mode which is asso- 
ciated with a single harmonic oscillator. According to 
quantum mechanics the energy of a harmonic oscillator 
with angular frequency lu is 

%) = M"+|) (4) 

where n is the excitation level. That is, the energy of the 
vibrational normal modes is quantized and n refers to the 
number of phonons in the crystal having the angular fre- 
quency to. Therefore, the following relationship between 
the excited classical amplitude spectrum of elastic waves 
and the phonon population holds 

Afe) 2 oc(n(^) + i)M«^l (5) 
z oji 

in which index i identifies the normal mode. The combi- 
nation of (2) and (5) yields 

(W)) t -MQ) /|^j (6) 

which is a one-dimensional equivalent of the relation for 
TDS derived by Warren 36 . 

We thus conclude that UXRD from a quasi- 
monochromatic strain pulse directly measures the 
squared spectral amplitudes of the plane elastic waves 
constituting the strain pulse. As an example we consider 
the reciprocal lattice vector G o2 of STO and rewrite (2) 
and (6) into 

A(q) oc 002 + q))t - iu P (Goo2 + q) (7) 

n(q) oc uj(q)((I p (G 00 2 + q))t - I up (G 02 + q)) (8) 

In the standard 9-29 geometry applied in our UXRD 
experiments, the magnitude of the phonon wave vector q 
is: 

4-7T 

q= — | sin 6» -sin 6» 1 (9) 
Ax 

where Ax is the X-ray wavelength, 9 is the X-ray in- 
cidence angle with respect to the sample surface ((001) 
crystallographic plane) and 9q is the Bragg angle. 
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FIG. 2. (Color online) Comparison of the calculations of the 
spectral strain amplitudes of phonons with the corresponding 
x-ray intensity difference signal of perturbated-unperturbated 
structures. The red solid line indicates the calculated phonon 
spectrum the sample in the q domain. The values in the verti- 
cal axis correspond to the phonon spectral amplitudes calcu- 
lated for the excitation fluence of 6 mj/cm 2 . The blue bullets 
show the square root of the x-ray intensity difference signal 
calculated for the pertur bated and unpertur bated structures. 
The vertical scale factor for the UXRD signal is arbitrary. 



To demonstrate, that (7) is applicable to our case we 
compare the calculated amplitude spectrum of the laser- 
excited strain waves to the dynamical UXRD simulations 
from the same acoustically perturbed sample 20 . The 
calculated spectrum for the laser fluence of 6 mJ/cm 2 
is plotted in Fig. 2 as a red solid line. The dynami- 
cal UXRD calculations are performed for 200 time steps 
within the interval from 100 ps to 300 ps after the exci- 
tation and then time-averaged. The blue symbols in the 
Fig. 2 show the scaled time averaged square root of the 
intensity differences [cf. (7)] obtained from the dynami- 
cal UXRD calculations. 

We see that the curves almost coincide although the 
fine structure of the UXRD-related curve (blue sym- 
bols) slightly deviates from the actual spectrum (red solid 
line). This is due to the fact that the scattered X-ray 
intensity from a propagating strain wavepacket actually 
oscillates in the time domain at each fixed q-vector with 
the frequency of the corresponding phonon mode. These 
oscillations were successfully observed in pioneering ex- 
periments with the advent of UXRD 33-35 . The classical 
explanation for the oscillations is the interference of X- 
rays scattered from a moving grating (elastic wave) and 
the X-rays scattered from the static component of the 
crystal lattice. This is an X-ray analog of Brillouin os- 
cillations in all-optical experiments 37 . In our case, as we 
can see, this interference is not very strong. Therefore, 
the shape of the blue curve slightly depends on the aver- 
aging time window. To eliminate these artifacts, the av- 
eraging time window should be either much longer than 
any phonon vibration period or we need to fit an inte- 



ger number of vibrations for each q. The averaging over 
many vibrational periods is not possible in our case, be- 
cause the actual phonon life time is only several vibration 
periods as we will see later. 

To finish this section we briefly review the conditions 
at which the approximation (7) should be valid: 

• The interatomic displacement in the strain wave is 
much less than the interatomic distance: 

\r m -r n \ < m7 in \a(m- n)\ (10) 

in which a is the interatomic distance, m and n are 
the index number of atoms. This is required by 
the perturbation theory of X-rays scattered from a 
dynamical lattice 36 . 

• The wave vector of a phonon is much smaller than 
any reciprocal lattice vector: 

q<G (11) 

This is necessary to avoid the signal overlap from 
the adjacent Bragg reflections of the crystal 

• The wavelength of a phonon mode is much smaller 
than both the X-ray coherence and the X-ray ex- 
tinction lengths: 

q » y (12) 

in which I is either X-ray coherence or X-ray ex- 
tinction lengths, depending on which one is larger. 

• The X-ray intensity is time-averaged over many vi- 
brational periods 

III. EXPERIMENTAL RESULTS 

We performed UXRD experiments on the laser-excited 
five-period SRO/STO multilayer with the structure pa- 
rameters presented in the Table I. In this section we 
present the experimental results which evidence the pres- 
ence of a propagating quasi-monochromatic LA phonon 
wavepacket. We discuss the dynamics of the first and 
second-order transient diffraction peaks and the corre- 
sponding dynamics of the strain pulse. 

In the experiment we acquire the X-ray photons scat- 
tered from the sample 50 ns before each optical pulse and 
at a given probe delay after each optical pulse. The cor- 
responding scattered X-ray intensities from perturbated 
and unperturbated sample are thus defined as I p and I up , 
respectively. The time resolution of the experiments was 
100 ps due to the limited X-ray pulse length, therefore the 
measured X-ray intensity is time averaged over multiple 
phonon vibrations. 

Figure 3(a) shows the measured x-ray intensity differ- 
ence signal I p — I up (blue symbols) in the vicinity of the 
STO (002) Bragg peak (not visible) at 100 ps after the 
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FIG. 3. (Color online) (a) The solid dots show the UXRD 
difference signal with an arbitrary vertical scale factor. Red 
is the calculated spectrum of the synthesized wavepacket. (b) 
The integrated UXRD difference signal of the first two order 
phonon peaks as a function of time. The solid lines indicate 
the exponential fit. 



laser-pulse excitation. The experimental incidence angle 
has been converted into the phonon wavevector using 
(9) . The vertical scale for the measured data is arbitrary. 
The UXRD data exhibit the first and second-order spec- 
tral components of the synthesized quasi- monochromatic 
phonon wavepacket inside the STO substrate at wave 
vectors q\ « 0.045 rad/nm and (72 ~ 0.09 rad/nm, respec- 
tively. Given the longitudinal sound velocity in STO (cf. 
Table I) , the linear phonon dispersion relation of acoustic 
phonons implies the corresponding hypersonic frequen- 
cies v\ 55 GHz and « 110 GHz. The non-vanishing 
contributions between the phonon peaks are due to the 
diffraction from the laser-heated multilayer. However, 
since the lattice constants throughout the multilayer are 
larger than that of the substrate (cf. Table I), the X- 
ray scattering from the multilayer is rather weak in this 
angular range. The red solid line in Fig. 3(a) shows the 
squared amplitude spectrum of the propagating sound 
wave as obtained from the linear-chain model. The shown 
spectrum includes the convolution with a Gaussian res- 
olution function having a full width at half maximum 
(FWHM) of 15 • lO^nm" 1 to fit the angular resolu- 
tion of the UXRD experiment. The main contribution 
to the XRD peak broadening is due to the sample bend- 
ing according to the stationary laser heat load 13 . The 
UXRD signal shows very good agreement with the con- 
voluted spectrum in terms of position, relative intensity 
and width of the first and second-order phonon peaks. 
This verifies the relation between the measured x-ray in- 
tensity and the amplitude spectrum of the coherent strain 
wave derived in (7). 

In the following we discuss the intensity changes of 
the measured phonon peaks with time. During the first 
100 ps after laser-excitation the intensity of the phonon 
peaks builds up 16 due to the unfolding of the initially ex- 
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oxp , the apparent decay time 



TABLE II. Comparison of 
UXRD intensity decay time t< 
due to x-ray absorption r a b s and the derived phonon lifetime 
r p h for the first and second-order phonon peaks. The corre- 
sponding standard deviations <7 cxp / abs are also shown. 
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0.49 ± 0.13 



cited superlattice phonon mode into the substrate 1,21 . 
Subsequently, the integrated intensity of the phonon 
peaks decays exponentially as is evidenced by the blue 
symbols in Fig. 3(b). The red solid lines show fits ac- 
cording to the function 



f(t) = AI e 



' cxp 



(13) 
arc the 



where the two fitting parameters AIq and r, 
amplitude and decay time of the measured signal. The 
data points before 200 ps after the excitation were ex- 
cluded from the fit since the wavepacket may not yet 
fully propagated from the multilayer to the substrate. 
The extracted decay times T exp and standard deviations 
<7 exp are shown in Table II. 

There are two major reasons for the observed decrease 
of the UXRD peak intensities. First, the absorption of 
the X-rays by the crystal reduce the sensitivity of the 
X-rays to the wavepacket as it propagates deeper into 
the substrate. Second, the dissipation of energy from the 
elastic wave due to the finite phonon lifetimes leads to a 
decay of the strain amplitude. 

Considering the first reason, the decay time of the 
UXRD signal exclusively due to X-ray absorption is re- 
lated to the X-ray absorption coefficient \i — 0.056 [imT 1 
by 
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(14) 



where c s = 7.9 nm/ps is the longitudinal sound speed in 
the substrate 26 . The relative variation of the x-ray scat- 
tering vector Q during the presented UXRD experiments 
is 10~ 3 which implies that r a b s is virtually independent 
of the observed phonon wavevectors q. In the measured 
off-Bragg region the X-ray extinction due to dynamical 
X-ray diffraction is negligible, therefore only the angular- 
independent X-ray absorption is relevant. Under the cho- 
sen experimental conditions we estimate a signal decay 
time of T a bs ~ 466 ps due to the x-ray absorption. Never- 
theless, since this value is critical for the correct interpre- 
tation of the experimental data, we have performed dy- 
namical XRD calculations based on results of the linear- 
chain lattice dynamics in harmonic approximation which 
excludes the effect of phonon damping 16,21 . The simula- 
tions yield the g- independent value of 450±5ps for the 
decay constant due to X-ray absorption which we will use 
in the following. 
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FIG. 4. (Color online) (a) Transient first and second-order 
phonon diffraction peaks at various probe delays. The in- 
tensity decay due to sound attenuation is accompanied by a 
continuous shift to larger wave vectors as is indicated by the 
arrows, (b) First and second-order peak position determined 
by Gaussian fits as a function of time delay. The solid lines 
show the linear fits to the experimental points. 



Regarding the second reason for the decay of the 
UXRD phonon signals, we assume an exponential law 
for the decrease of the phonon population n(g, t) and de- 
fine the associated decay time r p h- According to (7) the 
corresponding intensity of the scattered X-rays possesses 
the same decay constant. 

Therefore, the UXRD signal decay mechanisms intro- 
duce the relationship between the phonon lifetimes r p h, 
X-ray absorption time constant r a b s and the experimen- 
tally measured time constant r oxp : 
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(15) 



The phonon lifetimes r p h extracted from the UXRD 
experiments according to (15) are given in the Table (II) 
for the wave vector magnitudes corresponding to the first 
and second order phonon peaks. The standard deviation 
<Tph for the phonon lifetimes are calculated from the stan- 
dard deviations <7 xp of the experimental time constants 
according to the error propagation relation. 

The important observation is the fact that the deter- 
mined lifetimes differ by a factor of approx. 4 whereas 
the related phonon wave vector differ by a factor of 2. 
In addition to the signal decay we see a gradual drift of 
the phonon peak position to higher values of q as the 
wavepacket propagates. This is seen from the Fig. 4(a) 
in which the phonon spectrum is shown for different time 
delays. 

In the following section we discuss the physical inter- 
pretation of the described observations. 



IV. DISCUSSION 

The central observation from the transient phonon- 
induced diffraction peaks presented in the previous sec- 
tion is the apparent quadratic decrease of the phonon life- 
time with phonon wave vector q. This observation is in 
agreement with the 1/w 2 law predicted by the Akhiezer's 
sound attenuation mechanism 38 . 

There are basically two theories which explain the at- 
tenuation of a hypersonic waves in dielectric crystals 
due to incoherent anharmonic phonon-phonon scattering. 
These are the Landau-Rumer theory 39 and the Akhiezer 
theory 38 which have different application limits: 
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where uj is the angular frequency of the hypersonic wave 
and T t h is the mean thermal phonon relaxation time. 
That is, in both cases the elastic energy of hypersound 
waves decreases with time due to the interaction with 
incoherent thermal phonons. The mean thermal phonon 
relaxation time Tth can be estimated from the thermal 
conductivity k, heat capacity Cy, average sound speed 
v s and mass density p by the relation 27 



k = l^pCvvlnh 



(18) 



For STO at room temperature (18) yields T t h ~ 0.26 ps. 
The experiments are performed at room temperature. 
More precisely, the sample persisted at around 400 K 
during the actual measurements due to the thermal load 
from the laser 13 . Condition (16) is fulfilled for this tem- 
perature range. That is, for the presented experiments 
Akhiezer's theory of relaxation damping could be ap- 
plied, hereby explaining the observed ratio of the phonon 
lifetimes for the first and second-order phonon diffraction 
peaks. However the UXRD data exhibit additional fea- 
tures which cannot be explained within Akhiezer's sound 
attenuation model. We extracted the transient peak posi- 
tions by Gaussian fits to the data shown in Fig. 4(a) and 
the results are plotted as symbols in Fig. 4(b). The solid 
lines indicate linear fits to the phonon peak positions as 
a function of time. Clearly, a gradual shift of the first 
and second-order spectral components to higher g-valucs 
can be observed as the phonon wavepacket propagates 
deeper into the STO substrate. 

A recent study revealed the nonlinear propagation 
of large-amplitude sound wavepackets in STO at room 
temperature 23 . In addition to the Akhiezer- like attenu- 
ation of the coherent LA phonons the authors observed 
transient changes of the acoustic spectrum due to co- 
herent anharmonic phonon-phonon scattering within the 
wavepacket. The lattice anharmonicity gave rise to a 
strain-dependent longitudinal sound velocity. In partic- 
ular, the sound velocity of compressive (tensile) parts of 
the wavepacket was found to increase (decrease) with the 
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strain amplitude. This effect led to an anomalous disper- 
sion of the wavepacket and the corresponding modifica- 
tion the phonon spectrum. 

Accordingly, we expect the first compressive half of 
the wavepacket shown in Fig. 1(b) to propagate faster 
than the second tensile half. Moreover, the individual 
pulses inside the respective parts also exhibit different 
velocities due to the exponential amplitude distribution 
determined by the optical penetration of the pump light 
in SRO. For both the compressive and tensile parts of 
the strain pulse the spatial separation A of the individual 
pulses of the wavepacket is reduced as it propagates, i. 
c. the wavelength of both sub-packets is decreased. This 
explains the observed shift of the phonon peaks to larger 
q- values. 

The presented UXRD data thus evidence the influ- 
ence of two different effects on the propagation of LA 
phonon wavepackets generated by periodic multilayers. 
First, the inevitable attenuation of the wavepackets by 
Akhiezer's relaxation damping and, second, the change of 
spatial and spectral shape of the wavepacket by nonlinear 
sound propagation. Both effects influence the observed 
phonon lifetime, however, for a quantitative determina- 
tion of the respective contributions additional measure- 
ments have to be performed. The results of our earlier 
all-optical experiments having much stronger excitation 
were successfully explained solely in the framework of 
nonlinear acoustics 23 . We believe that at the presented 
experimental conditions the influences of both damping 
mechanisms the non-linear acoustic propagation and the 
Akhiezer's relaxation arc comparable. 

UXRD has the advantage of measuring the lattice dy- 
namics directly and quantitatively, i.e. the absolute am- 
plitude of the lattice motion is determined. The wavevec- 
tor range over which acoustic phonons in bulk material 
are accessible is very large. In particular resolving the 
second order phonon peak as presented in this paper or 
higher orders is possible. The extension of the UXRD 
detection of acoustic phonons in amorphous materials is 
challenging the available x-ray fluence, since the Bragg 
spots are dispersed in diffraction rings. 

On the other hand all-optical picosecond 
acoustics 3 ' 4 ' 7 ' 8 ' 17 , in principle, do not require crystalline 
materials and for transparent media, the propagation 
of strain pulses can be monitored over longer distances. 
With current technology femtosecond time resolution 



is standard in all-optical experiments, while it is still a 
challenge in x-ray technology, which essentially resolved 
by free-electron lasers. High time-resolution permits 
the determination of the wave vector dependent sound 
velocity in addition to the damping time. 

We believe that the UXRD-bascd methods and the 
all-optical methods do not compete with each other but 
complete each other, together providing a more complete 
picture of the complex coherent phonon dynamics for a 
broader range of frequencies and wavevectors and for a 
broader class of materials and experimental conditions. 
V. CONCLUSIONS 

This report presents ultrafast x-ray diffrac- 
tion (UXRD) studies on laser-excited periodic 
SrRuOs/SrTiOa multilayers which are epitaxially 
grown on a SrTiC>3 substrate. The ultrafast heating of 
the metallic SrRu03 layers by ultrashort laser pulses 
generates coherent longitudinal acoustic phonons which 
eventually propagate into the substrate as a quasi- 
monochromatic coherent LA phonon wavepacket at 
hypersonic frequencies. We discussed the properties 
of such wavepackets in detail and derive equations 
which show that UXRD is a powerful tool to measure 
the spectral phonon population and its dynamics. 
The presented UXRD data evidence the formation of 
quasi-monochromatic coherent phonon wavepacket. We 
extract the phonon lifetimes of the first and second-order 
peaks of the phonon spectrum. The observed quadratic 
decrease of the phonon lifetime with increasing phonon 
wavevector q is in accordance with Akhiezer's mechanism 
of relaxation damping. Shifts of the peaks corresponding 
to the excited phonons to larger q-values are interpreted 
as a modification of the spatial shape profile due to 
the non-linear wave propagation leading to a strain 
dependent sound velocity. This considerably modifies 
the observed phonon lifetimes. 

In essence, UXRD provides a detailed and direct 
view on the complex nonlinear evolution of phonon- 
wavepackets, including incoherent damping of the 
phonon amplitude by coupling to other modes and spe- 
cific coherent changes of the wavevector spectrum. 
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